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q^ \ We study out of equilibrium thermal field theories with switching on the interaction 

<3\ ' occurring at finite time using the Wigner transforms of two-point functions. 

For two-point functions we define the concept of projected function: it is zero if any of 
times refers to the time before switching on the interaction, otherwise it depends only on 
Q_ij the relative coordinates. This definition includes bare propagators, one-loop self-energies, 

etc. For the infinite-average-time limit of the Wigner transforms of projected functions we 
define the analyticity assumptions: (1) The function of energy is analytic above (below) 
the real axis. (2) The function goes to zero as the absolute value of energy approaches 
/\ ■ infinity in the upper (lower) semiplane. 

Without use of the gradient expansion, we obtain the convolution product of projected 
functions. We sum the Schwinger-Dyson series in closed form. In the calculation of 
the Keldysh component (both, resummed and single self-energy insertion approximation) 
contributions appear which are not the Fourier transforms of projected functions, signaling 
the limitations of the method. 

In the Feynman diagrams there is no explicit energy conservation at vertices, there is 
an overall energy-smearing factor taking care of the uncertainty relations. 

The relation between the theories with the Keldysh time path and with the finite time 
path enables one to rederive the results, such as the cancellation of pinching, collinear, 
and infrared singularities, hard thermal loop resummation, etc. 
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1 Introduction 

Out of equilibrium thermal field theory [I], |2j has recently attracted considerable interest |3], [|, 

In many applications one considers the properties of the Green functions of almost equili- 
brated systems, at infinite time after switching on the interaction. A recent approach based on 
first principles has been successful in demonstrating the cancellation of collinear Jl6|, [HJ and 
pinching singularities [18, [19], GO, ETl B2, 123], the extension of the hard thermal loop (HTL) 



3 m m, E 

approximation |24|, |2^, |2(| to out of equilibrium |27, 28| , and applications to heavy-ion collisions 
29], RD], plj. A weak point of the approach was that most of these results were obtained under 



the assumption that the variation of slow Wigner variables could be ignored or, in other words, 
that these results were valid only in the lowest order of the gradient expansion [EJ2, 53, [M . 

For some problems, e.g., heavy- ion collisions, both limitations are undesired. One would 
like to consider large deviations from equilibrium. One cannot wait infinitely long as these 
systems go apart after a very short time, probably without reaching the stage of equilib- 
rium. In nuclear collisions, short-time-scale features have been studied in a number of papers 
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In the theories with switching on the interaction infinitely long before the time of interest, 
one tries to get some information by extrapolation to early times. However, in doing so, this 
information is either deformed or lost. Indeed, relaxation phenomena include many processes 
that are expected to terminate as one goes sufficiently far from the switching-on time of the 
interaction. Thus one can expect that the full theory also describes early times. In many 
papers, problems with finite switching-on timet!, especially those related to the inflatory phase 



of the Universe (see Ref . 1 1 1 ] for further references) have been studied by Wigner-transforming 
in relative distance and studying the dependence on relative time directlyCrl. In such an 
approach, one relies heavily on differential equations and numerical methods. 

In an attempt to remove weak points in both cases, we suggest the application of the method 
of Wigner transform (but now also in relative time) to the case of switching on the interaction 
at finite time. The time-integration path C (see Fig. 1) is now closing the part of the real axis 
between £; (the switching-on time) and tf (the switching-off time); in the rest of the paper the 
switching-off time is pushed to infinity. If we push the switching-on time to minus infinity, the 
connection to the Keldysh integration path is established. It turns out that this connection is 
highly nontrivial. 

To understand the limitations coming from the finite switching-on time, we start with a 
generic two-point function G(x,y). The quantities x and y are four- vector variables with the 
time components in the range tj < xo,yo < oo (here ti is the time at which we switch on the 
interaction; it is usually set to — oo, but we set it to £j = 0). We define the Wigner variables s 
(relative space-time, relative variable) and X (average space-time, slow variable) as usual: 

v x + y 

X = — — , s = x-y, 



G(x,y) = G(X + ^X-^ 



The lower limit on xq, y^ implies the following conditions on X Q and Sq. < X , 
s < 2X . To define the projected function (truncated, "mutilated function 
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(1.1) 

2A < 
PF in further 



text), we add two additional properties: the function does not change with X (homogeneity 
assumption), it is a function of (sq, s) within the interval — 2X < s < 2X Q and is identical to 
zero outside the interval: 

F(X + |, X - |) = e(X )6(2X - s )Q(2X + s )F(s , s), 

F(a ,3)= lim F(X + t X -±). (1.2) 

Xo— >oo Z Z 

Analogously, the Wigner transform of the projected function (WTPF) is obtained from the 
Wigner transform of the function defined on the infinite carrier of Sq (—00 < Sq < 00) with the 
help of the projection operators Px , which are the Fourier transforms of the above given 0's: 



/oo 
d P0 P X (PO, Po) F oo (p' ,p), 
-00 



Px {Po,Po) = -@(^o) ; , hm P Xo {Po,Po) = S{Po-Po), (1-3) 

^ PO ~ Po Aq^oo 



where the subscript " 00" , as it appears in Foo in (|1.3|) and in other expressions of this paper, 
is the short notation for the "limx ^oo"- As an illustration, the Wigner transform of the 
convolution product of two-point functions 

C = A * B & C(x, y)= f dzA(x } z)B(z, y) (1.4) 

is given by the gradient expansion (note that we have assumed the homogeneity in space 
coordinates, which excludes any dependence on X): 

Cx (po,p) = e-*<>A Xo (p ,p)B Xo (p ,p), = \{d Xo d* - d*d% ). (1.5) 

This general expression should be contrasted with our result valid for A and B being projected 
functions (see Sec. 5): 

r P01 + P02 1 ie~ iX °( p01 ~ P02+i ^ 

Cx (po,p)= dp 01 dp 02 Px (,Po, = )t: -^A^pox.pjB^p^p). (1.6) 

J 2 2tt p 01 - p 02 + te 

The analytic properties of the WTPF in the Xo — ► 00 limit as a function of complex energy 
are very important for further analysis. We define the following properties: (1) the function 
of po is analytic above (below) the real axis, (2) the function goes to zero as |po| approaches 
infinity in the upper (lower) semiplane. The choice above (below) and upper (lower) refers to R 
(A) components (note here that we do not require such properties for the Keldysh components). 
It is easy to recognize that the properties (1) and (2) are just the definition of the retarded 
(advanced) function. However, it is nontrivial, and not always true, that the functions with the 
R (A) index satisfy them. 

Under the assumption that A or B satisfies (1) and (2) (A as advanced or B as retarded) 
Eq. ( |1.6| ) can be integrated even further. We obtain 

C Xo (Po,P) = J dp' P Xo (po,Po) A oo(po,P) B oo{Po,P)- (1-7) 



The convolution product of two two-point functions which are WTPF's and satisfy (1) and (2) 
is also an WTPF. This product is then expressed through the projection operator acting on a 
simple product of two WTPF's given at Xq = Xq = oo. In the X — > oo limit the convolution 
product of two WTPF's which satisfy (1) and (2) is equal to the lowest-order contribution in the 
gradient expansion. The result makes sense since moving Xq — > oo is equivalent to ti — > — oo 
at fixed Xo in the standard approach. 

We find that some quantities, obtained at low orders in the perturbative expansion, e.g., bare 
propagators, one-loop self-energies, belong to WTPF's. This enables us to sum the Schwinger- 
Dyson series with the propagators and self-energies as both these quantities are WTPF. Under 
the conditions (1) and (2), the retarded, advanced, and Keldysh components at finite Xo 
are obtained by a simple action (smearing) of the projection operator onto the corresponding 
quantities obtained at Xq = oo, and the convolution product is a simple multiplication: 



Gr,x (Po,P) = J dpoiPx (pa,Poi)GRoo(poi,P), 

and similarly for the advanced component. The calculation of the Keldysh component requires 
a more elaborate treatment: one reduces the multiplication to a double (a single in some cases) 
convolution product; the result contains terms that are non-WTPF terms. Also, the one-loop 
contribution to the Keldysh component is corrected by non-WTPF terms. 

Now, it is Gr,oo(po,p) as given in Eq. ( |1.8| ) (and similar expressions for the advanced 
and Keldysh components and the single self-energy insertion approximation to Gr(a,K),<x>) to 
which previous results on the cancellation of pinching singularitiesliSr Ii3 Eii c3 and the HTL 



resummat ion ell' c-3 (and also on the cancellation of collineartLHl LL.LJ and infrared singularities 
if the properties (1) and (2) hold at the two-loop level) apply. 

From our study one can deduce a general rearrangement of the perturbation expansion at 
the non- Keldysh integration path: the contributions look like the zeroth order of the gradient 
expansion with the slow coordinate (Xo) pushed to +oo, but the use of the PF manifests itself 
as the appearance of the (J2jQoj + ^e)" 1 factor instead of — ttS(J2j Qoj) at each vertex, and as 
an overall projection (smearing) operator instead of the exact conservation of energy. 

Our study suggests that the results obtained by using the Keldysh integration path (ti — > 
— oo) could be related to the results of our approach (ti finite). This relation is possible at low 
orders of the perturbation expansion, i.e., as long as the expressions involved are the projected 
functions not breaking assumptions (1) and (2). Technically, the amplitudes are related by 
( p.. 3D , where " Fx " is the contribution in our approach (ti finite) and "Fqo" is substituted 
by the corresponding lowest-order contribution in the gradient expansion in the theories with 
ti — ► — oo. 

Finally, we note that our method is helpful in problems related to the time evolution of the 
system. Additional problems related to the gradient expansion in space components, appearing 
together with space inhomogeneity, will not benefit from our method. 

The paper is organized as follows. In Sec. 2 we give a general setup of out of equilibrium 
thermal field theory. In Sec. 3 we define finite-time Wigner transforms, define the projection 
operators, and introduce the notion of projected functions. We define analyticity assumptions 



(1) and (2). In Sec. 4 we define a few important examples of projected functions namely bare 
propagators and one-loop self-energies, and find that they satisfy the analyticity assumptions 
(1) and (2). In Sec. 5 we analyze the properties of the product of two and n two-point 
functions. We discuss the reduction of the inverse bare propagator to the space of projected 
functions. In Sec. 6 these properties are used to study the product of two pole contributions, 
and to discuss the reduction of the inverse bare propagator to the space of projected functions 
and the equations of motion. We sum the Schwinger-Dyson series for retarded, advanced, and 
Keldysh components of the propagator. We discuss the appearance of pinching singularities 
in our scheme. Section 7 is devoted to the modifications of Feynman rules in coordinate and 
momentum space. It is indicated that, in the absence of breakdown of assumptions (1) and 
(2), all energy denominators appearing at vertices can be replaced by delta functions. Section 
8 is a summary of the results and ideas described in the paper. 

2 Out of Equilibrium Setup 

We start by assuming that the system has been prepared at some initial time U = (To avoid 
inessential complications, we assume that the zero-temperature renormalization has already 
been performed). At ti the interaction is switched on and at time tf it is switched off (we shall 
take the limit tf — > oo). For ti < t\,t 2 < tf, the system evolves under the evolution operator 

I 

U{t 2 ,h) = T c [exp^•|dV£/(a; , )], (2.1) 

where c is the integration contour connecting t± and t 2 in the complex time plane and T c is the 
contour ordering operator. We provide T c with an extra property: for all times not belonging 
to the contour it gives zero. 

The Heisenberg field operator 4>(x) is obtained from the free field 4>i n (x) in the interaction 
picture as 

<f>{x) = U{ti,t)<j)i n {x)U{t,t i ), (2.2) 



(x 



T G [(j>in(x)ex.pi d 4 x'Ci(x')], (2.3) 

■J (.. 



where all fields on the right-hand side are in the interaction picture, and C is the integration 
contour of Fig. 1 (with the switching-off time pushed to infinity, tf —>■ +oo). In the Heisenberg 
picture, the average values of the operators are obtained as 

< 0(t) >= TrpO{t), (2.4) 

where p is the density operator admitting the Wick decomposition. Specially, we define the 
two-point Green function as 

G {c \x,x) = -i< T c <f>(x)<f>(x f ) > . (2.5) 



With the help of ( [2.1|) it can be written as 

G (c \x,x r ) = -i<T c (exp[i^dV£j(a/)]0 in (x)<M^ > • (2-6) 

In Eq. ( p.6|) it is implicitly assumed that the interaction Lagrangian does not depend on time 
explicitly. Indeed, direct time dependence through time-dependent perturbation, or through 
the background field which depends directly on time, ip(x, t) = (p{t) + ip(x, t), would break our 
scheme through the appearance of two-point functions which are not projected functions. 

We assume the single-particle density operator to be stationary with respect to the free 
Hamiltonian Hq = J2j Hj (f° r an alternative choice of the initial density, see Ref. 



P = X) IV>n > Pn < -0n| = -= exp(-^^/f i ), (2.7) 

Z j 

where the "temperature" function tl (3j (the "temperature" of the j th degree of freedom) is 
adopted to obtain the given initial-state particle distribution. Now, one has 

]T PjHj = J d A p[3( Po ) Po e( Po )5( P 2 -f- m 2 )a;a p , (2.8) 

j 

and obtains the distribution function f(po) as a function of (3(po), 

1 1 

UiPo) = —Trp p = —— — — -, (2.9) 

Po exp/?(p )poTl 



or the inverse relation 



log(-^ ± 1) 

PM = -^^ — -■ (2-10) 



Po 
The free fields are expanded in creation and annihilation operators as 



^-Swk; (a ' e ~ ,r ' +a ' e ' n ' (2 - n) 

with p = uj p = (p 2 + m 2 ) 1 / 2 : 

< a+a p > >= (2n) 3 2uj p f (uj p )5(p - p) , 

< a p a+ >= (2tt) 3 2^(1 ± f(u p ))6(p- j?), (2.12) 

where f(oo p ) is the given initial distribution. 

For completeness, in Eqs. fl2.8|) -( pT2 ) we include free fermions (lower case) in addition to 



free bosons (upper case). In the case of spin-1/2, spin-1, or higher-spin particles, additional 
spinor or tensor indices must appear. For simplicity, we do not show them explicitly. 
The noninteracting contour Green function is given as 

Gf n \x,x') =-%< T c in (x)0 in (x') > . (2.13) 

Depending on whether the times xq and x' belong to the upper ("1") or lower ("2") part of the 
path C, the function G\ n (x, x') splits into the components G^ Vi i n (x,x), //, v = 1,2. For 
the times xq < or x' < 0, the Green function is equal to zero owing to our definition of Tq. 



3 Projected Functions 

Let us start with the two-point function G(x, y). The quantities x and y are four- vector variables 
with time components in the range tj < Xq, y^ < oo (here tj is the time at which we switch on 
the interaction; it is usually set to — oo, but we set it to ti — 0). We define the Wigner variables 
as usual: 

v x + y 

X = ~^--, s = x-y, 

G(x,y) = G(X+ S -,X- S -). (3.1) 

The lower limit on xq, yo implies < Xq, — 2Xq < sq < 2Xq. The values of the function G 
for the (X, s) not satisfying these conditions are physically irrelevant. Our definiton of time 
ordering operator (see Sec II) sets them to zero, so that we can rewrite Eq. (|3.1|) as 



G(x, y) = 0(X O )0(2X O - s )Q(2X + s )G(X + ~, X - |). (3.2) 



Note here that the function G defined by Eq. (|3.2|) in general depends on Xq. At the points 
(X, s) which do not belong to the carrier of projection operator, the values of G are arbitrary. 
This freedom is used to define projected functions. The two-point function can be expressed in 
terms of the Wigner transform (i.e. Fourier transform with respect to Sq, Sj): 

G(X + ^,X- S -) = (27r)- 4 |dV- l(poS °-^G(po,p;X). (3.3) 

Here 

dsQ / dhe^ s «-^G(X + -, X - -) 

-2X J 11 

= jH | ds J d 3 S e^ So -^e(Xo)0(2X o - sq)Q(2X + s )G(X + i,X - ~). (3.4) 

We adopt a simplifying assumption of the homogeneity in space coordinates. This assumption 
excludes any dependence on X and we drop it as an argument of the function. 

The product of functions is a projection operator with a simple Fourier transform 



2tc J-2XH n Pn-Vn 



1 — ^ f 2X ° dsoe^M = l e(Xo) sin (2X0^-^)) 

-2X 71 PO-Po 

and 



e- isoP 'oe(X Q )e(2X + sq)Q(2X - s ) = JdpQe- lSoPo P Xo (po,p ). (3.6) 

It is important to note that 

v p 1 >\ v l sin(2X (p -^)) 

km Px {Po,Po) = hm } = S(p -p Q ). (3.7) 

Xo-»oo Xq^oo n Po ~ Po 



There is a hierarchy of the Px projectors: 

Px om (Po,p"o) = J rfPo-Px (po,Po)-Px^(Po,p"o), X 0tM = min(X ,X' ). (3.8) 

In this paper, the projected function is a very special two-point function F(x,y) = F(X + 
s/2,X — s/2): it does not depend on X, it is a function of (s , s) within the interval — 2X < 
Sq < 2X and identical to zero outside: 

F(X + |, X - |) = 0(X o )6(2Xo - s )e(2X + s )F(s , s). (3.9) 

Function F is related to the limit X — > 00: 

lim F(X + |,X-|)= J P(s ,5). (3.10) 

Important property of the projected function is that the whole X dependence is introduced 
by the projection operator. 

/oo 
dp'oPx (po,Po)Foo(Po,P), (3.11) 

-00 

Important examples of projected functions are retarded, advanced, and Keldysh components 
of free propagators. Further examples will emerge in the next sections. 

For further analysis, the analytic properties of the Xo — >■ 00 limit of the WTPF as a function 
of complex energy are very important. We define the following properties: (1) the function of po 
is analytic above (below) the real axis, (2) the function goes to zero as \p \ approaches infinity 
in the upper (lower) semiplane. The choice above (below) and upper (lower) refers to R(A) 
components. 

One should note that the properties of the projected functions are tightly related to the 
abrupt cutoff at £». Any smoothing of the cutoff would also change these properties. 

4 Examples of Projected Functions 

4.1 Poles in the Energy Plane 

We start with the simplest projected functions: simple poles in the energy plane. The pole 
contribution to the Green function is 

Qoo,pole(Po) = — • (4.1) 

Po-Po 

For Impo < 0, it satisfies assumptions (1) and (2) as a retarded component, but not as an 
advanced component (and for Imp Q > 0, just the opposite). 
It can be projected to finite X : 

I _ e -2iX (po-po)sgn(Impo) 

Qx , P oie{po) = a = . (4.2) 

Po-Po 



For any finite X , the function Gx , P oie(Po) is regular at p . For large X , Eq. ( fO|) exhibits the 
exponential decay e - 2X o\ Im Po\ independently of the sign of Imp . 
It can be transformed back to the variables (X , s ): 

G P oie(X + y,X - y) = iae- lp0S0 (Q(s sgn(Imp )) - &(-s sgn(Imp ) - 2X )) . (4.3) 

Evidently, this contribution is a projected function. For Impo < 0, it is different from zero only 
at < so < 2X , i.e., it is a retarded function, and for Impo > 0, it is an advanced function. 

4.2 Propagator 



We start with Eqs. ( 2.11 ), ( 2.12| ), and (|2.13|) . The transition to the R/A basis is straightforward. 



Careful calculation gives for the retarded component (0 < xo, < y$): 

G,(., y) = -G,, + a v = / ^ p2 _ m 7 +2itpo e-^>, (4.4) 

and for the Keldysh component: 

G*(a;,3/) = Gi (1 + G 2 , 2 = J d A p2n5(p 2 - m 2 )(l ± 2/(w p ))e-^ x -^. (4.5) 

As our Gr and Gj<- depend only on s = x — y and vanish at times before switching on the 
interaction, they are pr ojected functions. The Wigner transform over the infinite XQ—yo interval 



gives as usualS \±3 LL^I (note, however, that we avoid^ using the nonanalytic function e(po) 
in the expression for Gk,oo)'- 

— ^ 
Gr,oc(p) 



p 2 — m 2 + 2iepo ' 

GkAp) = -(1 ± 2/K))^- 1 (p G R ^(p) - p G Aoo (p)) . (4.6) 

The e parameter, which regulates these expressions, should be kept uniformly finite during 
the calculations, and the limit e — > should be taken last of allQ. This specially means that 
limx ^oo ex p( — X e) = and the terms containing this factor vanish in the X — > oo limit. 
The finite Wigner transform (x > 0, y Q > 0, X Q = (x + y )/2) is obtained by smearing 

G R ,xM = [Px G RtOO ](p) = -G% Xo (p), G K , Xo (p) = [Px G K>oo ](p). (4.7) 

It is easy to verify that neither the spinor nor the tensor factor changes our conclusion (|4.7|) . 
One can even integrate expression ( f4.7|) . For a scalar particle, one obtains 

G°r,x (p) = 2 =g ~V - ( l ~ (cos2X uj p -z^sm2X uj p )e 2iX °^)) 
Po — P — m + 2«ep \ u p J 

= G°r,oo(Po) f 1 - (cos2X o; p - i^sin2X uj p )e 2tXo ^ +i A . (4.8) 

It is important to observe that, at any finite X , the above expression is not singular at p = 

±u p . 



Evidently, for X — > oo, the first term in Gr >Xo (p) gi yes Gr,oo, while the other two "oscillate 
out". For the Keldysh component, one needs 



.u n 



(poG R )x Q (p)=PoG% 00 ( P o) f 1 - (cos2X ^-^sin2X ^)e 2 ^ 0+ -M . (4.9) 

Then one can use the analogy to (|4.6|) : 

G°k,x (p) = "(I ± 2/K))^ 1 ((p G° R ) Xo (p) - (PoG° A ) Xo (p)) ■ (4.10) 

For a spinor particle, one obtains 

G%M = *7°e 24Xo(po+4e) K - ^) sin2X c p + G% Xo (p)(^p, + m). (4.11) 

Similarly, for a vector particle (for simplicity, we choose the Feynman gauge): 

GI u>R)Xo (p) = 9wG% Xo (p). (4.12) 



We note here that the explicit expressions (|4.8|) -( f4~T2J ) will not be necessary for further discus- 
sion. 

4.3 One-Loop Self-Energy 

To discuss the amputated one-loop self-energy, we start with (the underlying theory includes 
bosons and fermions with three point vertices, but spin and internal symmetry indices are 
suppressed for simplicity of presentation) 

£(x,y) = g 2 S(x,y)D(x,y) 

= ^|dV- ips ^ (po,^)5 00 (p[ ) ,^|Ae-^Px (go,^) J Doo(g( ) ,g1. (4.13) 

The Wigner transform (with respect to s = x — y) is 

£x„(Poi,Pi) = 9 2 r° ds J dp d 3 pdq e- l{po+qo)so e(2X - s )Q(2X + s ) 

xS oo {p' ,p)D oo (q' ,p 1 -P) 

= 9 2 J dp' d 3 pdq' P Xo (p 01 ,p' + q^S^ipo^D^q^p! -p), (4.14) 

where as an intermediary step we have used the representation of the bare propagators ( |4.7| ) 
and the representation of the projectors ( |3.5| ) and ( |3.6| ). Finally, one reads ( [4.14J ) in the R/A 
basis (S R(j4) = -(En + £21(12)), ^k = -(Sii + S 22 )), as 

Zr(A),x (p) = [Px Zr(A),co](p), ^k,x (p) = [P Xo ^k,oo](p)- (4.15) 
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To calculate £,r,oo and S^ j00 we start with Eqs. (11.23-25) in |23|1 . After taking into account 
that the product of the retarded with the advanced function, with the same time variables, 
vanishes, one obtains 

in f d \\ 
^R,oo(q) = -„- / j7r^j( h (ko, w fc ) + h(q - k Q , u q -k))D R>00 (k)S R>00 (q - k)F, (4.16) 

where D and S are bare scalar propagators, h(k ,ujk) = —koU> k (l ± 2f(u p )), and the factor 
F = F(k , \k\, q , \q\, kq, ...) includes the information about spin and internal degrees of freedom 
(F — 1 if all particles are scalars). Now, to verify that the one- loop self-energy E# )00 satisfies 
assumptions (1) and (2), one observes that the vacuum contribution satisfies them (for excep- 
tions, see, e.g., Refs. |£| ^3|), while the contributions to £_r i00 from various k points are linear 
and additive in distribution functions. 

For finite e, this contribution possesses singularities only below the real axis in the complex 
qo plane, and vanishes as |go| ~~ > oo in the upper semiplane. However, there is no guarantee 
that the imaginary part of Er i00 is negative. 

Using the same method one calculates the Keldysh component. Although assumptions (1) 
and (2) are not imposed on E^, one can decompose this component into two pieces T, K = 

Ex,i?(A),oo(g) = =F— J TTT^lt 1 + h(k ,u} k )h(q - k ,u q _ k ))D R(A)tOO (k)S R ( A)>oo (q - fc)F,(4.17) 

where ^k,r(A) satisfies assumptions (1) and (2) in the way the retarded (advanced) function 
does. General analytic properties of the expressions of the type ( 4.16 ) and (|4.17|) are well known: 



there is a discontinuity (cut) along the real axis, starting at thresholds for real processes and 
extending to ±oo. 

5 Convolution Product of Two Two- Point Functions 

Let us now consider the convolution product of two Green functions: 

C = A*B^ C(x, y)= f dzA(x, z)B(z, y). (5.1) 

In terms of Wigner transforms: 

C( Po ,p; X) = [ ° ds fd 3 s [d 4 ze i{poS °-P s) 

J-2X J J 



jL J dW-^^-^B^pi; A 2 ), 



(2tt) 
X 1 = X + ^, X 2 = X - ^, Sl = x - z, s 2 = z - y. (5.2) 
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The assumed translational invariance helps us easily integrate the space components of momenta 
and coordinates. To do so, we substitute d 3 sd 3 z by d 3 Sid 3 s 2 (Jacobian J = 1) 

s = Si + s 2 , z — h a . (b.d) 

The momenta should be equal (p = Pi = P2) and one obtains (note that the dependence on 
X's is reduced to the dependence on X 's; further in the text it is indicated as an index), 

Cx (po , P) = 77^ J ds I dz I dp 01 I dp 02 

xe^^-^^-^^Ax^ipox^Bx^ip^p). (5.4) 

For energy integrals, we proceed in a somewhat different way. We shrink our choice of the 
functions A(x, z) and B(z, y) to the projected functions. Then we can use the connection to 
the Wigner transforms on the infinite carrier: 

C Xo (Po,P) = ^/_ 2 2 ^ rfs o/^o/rfpoi/^o 2 e i(poSo - poiSoi - po2S02) 

x J dp' ol Px 01 ij)'Qi 1 Pox)A 00 {pQ 1 ,p) J dp' 02 Px 02 (po2,Po2)B 00 {p 2,p). (5.5) 

The integration dp' 0l dp' 02 is easily performed with the help of Eq. (|3.6|) and one obtains 
C Xo (po,p) = 7^2/^ ds J dz J dp 01 J dp Q 2e^ s °- po ^- p ^ s ^ 

x6(2X 01 + soi)e(2X i - s i)e(2X 02 + s 02 )e(2Xo2 - s 2)A oo (p 01 ,p)B oo {p 2,p). (5.6) 
The product of 6 functions is transformed into 0(2X O + so)6(2X — so)O(zo)- Then 

Cx (Po,P) = J J dpo 1 dp 02 8(po,poi,Po2)A 00 (p 01 ,p)B 00 (p 2,p). (5.7) 



Here 



5(p ,Poi,Po2) = 7^ r° ds * i" ^ e i(poSo ~ PoiSol ^ 02S02) 
(27r)^ J-2X Jo 



1 r2X 



/ZJi roo 

dso / dzoe 
-2Xn JO 



= i( so (p _£Qi+£02) + ( Zo -Xo)(p 01 -po2+ie)) 



(27r) 2 J-2X Jo 

= Px (po, ^±^)^ •__ e -«b(« tt -*»^), (5.8) 

2 2tt poi - P02 + te 

where we have used 

/ dz e lza = . (5.9) 

jo a + le 

We can write the final expression as 

r P01+P02 1 ie~ iX ° < - P01 ~~ P02+ie ' ) 

Cx (po,p)= dpoidp 2Px (po, = )» -^A 00 (p 01 ,p)B 00 (po2,p}- (5.10) 

J 2 2tt poi - P02 + ie 
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Expression ( p. 10 ) is the key for finite-time thermal field theory. 



If A is an operator satisfying assumptions (1) and (2) for advanced components, we can 
integrate expression (|5.10|) even further. After closing the poi integration contour in the lower 
semiplane, one obtains (if B is an operator satisfying (1) and (2) for retarded components, one 
can achieve the same result by closing the P02 integration contour in the upper semiplane): 

Cx {po,P) = J dp iPx {po,Poi)A OD {poi,p)B 00 (p 01 ,p). (5.11) 

This is an extraordinary result: the convolution product of two WTPF's is a WTPF under 
conditions (1) and (2). 

As expected, in the X Q = oo limit, Eq.( |5.1lD becomes a simple product 



lim C Xo (po,p) = Aooipo^BooipoiP). (5.12) 

Xq— »oo 

At finite X , Eq.( |5.11|) exhibits a smearing of energy (as much as it is necessary to preserve 
the uncertainty relations). 

5.1 Convolution Product of n Projected Functions 

The product of n two-point functions is obtained by repeating the above procedure: 

~ n—l 

Cx {Po,P) = J n(dpa,j)dpo, n Px (Po, (Po,i+Po,n)/2) 

x n (^i,oo(Po,i,j5)^- -A e- iXo(p ^- po '^ n - 1)e) A nt00 (p 0jn ,f). (5.13) 



For the intermediate products in Eq. ( 5.1 3|) to hold we must require that at least n-1 of the 



functions in the product satisfy assumptions (1) and (2). Furthermore, the order of these 
functions is important: the retarded functions should be on the right, the advanced on the left, 
and the function that is neither advanced nor retarded in the middle. However, this is not the 
order in which the components appear in the Schwinger- Dyson equation. 

If the above requirement is fulfilled, one obtains (index R for the retarded component, a 
similar expression for the advanced component) 

1 % _ 1 i 



/i i 
d P0,j 7, -— A R,j,oo (po,j , P) 
27T On ,-_i —Vm+ie 



2vr poj-i - Po,j + ^ ' ' ' 27r Po,j ~ Poj+i + it 

1 i 

= A RJtOD (p 0d _ 1 ,p)- —. (5.14) 

27rp ,i-i ~Po,j+i + i£ 

Then one finds 

~ n 

Cx {po,P) = J dp OA P Xo {p ,Po,i) II A j,oo(Po,x,P)- (5.15) 
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6 Examples of Convolution Products 

6.1 Convolution Products of Pole Contributions 



We assume two pole contributions as shown in ( [4.1]) : A 00tP( j ie>i = a,i/(po — po,i), i — 1,2. The 



product C = A\ * A 2 is simple in the cases in which both contributions are retarded functions, 
or both are advanced, or A\ is an advanced and A 2 a retarded function. Then one can simply 
use ( |5.11| ) to obtain 



Cx (po) = dpoiPx (Po,Poi) — 2 ~^- (6-1) 

J PO- P0,1 PO - P0,2 

The case in which A\ is a retarded and A 2 an advanced function (i.e., Impoi < o, Imp 02 > 0) 
requires additional care. After substituting them into (|5.10|) , we choose new variables Pq = 
(Poi + Po2)/2 and A = p 01 — p 02 , and integrate over A to obtain 



Cxoipo) = / dP P Xo (po,Po)-^- 

J /o 



a\ a 2 

~ P0,1 Po ~ P0,2 



r I dP P Xo ( Po ,Po h ' op i ofp —^ ■ + Ofp - UT • ( 6 - 2 ) 

Poi + P02 ~ 2P \ 2 ( P o - Poi) - ie 2(P - P02) + *e / 



The first term is formally identical to ( |6.1|) . The second term consists of two non-WTPF pieces. 
As we shall see in the subsection discussing pinching singularities (Sec. 6.4), non-WTPF 
contributions appear also in the convolution product of the type Gr *Q* Ga or, more generally, 
in the convolution products containing retarded components positioned on the left from the 
advanced components (...*Ar*...*Ba*---)- The non-WTPF terms depend directly on X , they 
are carrying the nontrivial information about the time evolution (i.e., about the dependence 
on Xq). However, the non-WTPF terms cannot be convoluted further using Eq. ( [5.1 0| ) or 
( |5.11| ). This fact indicates the natural limits of the applicability of the methods developed in 
this paper. 

6.2 Inverse Propagator and the Equations of Motion 

To define the inverse propagator, we use the results of subsection (4.2). We define the restriction 
of Gft 1 on the subspace of projected functions as 

G RXo (PoiP) = J dp'oPxoipo^^G^ip'^p), G R ^(p ,p) = i(p 2 - m 2 + 2iep ). (6.3) 

This integral does not converge in the absolute sense, thus we cannot calculate the dependence 
of Gr 1 on Xq. Nevertheless, we can apply it from the left to some class of functions. For 
example, we can apply it formally to G\r,x : C/? 1 * Gr = 1, or written out more explicitly 

/ dp' P Xo {po,p )i{p 2 - m 2 + 2ip e)G R)OO {p' ,p) = 1. (6.4) 

This equality is obtained using a simple integration over p 02 in the expression of the type 
( p,10| ). We cannot verify the second identity Gr * G^ 1 = 1 directly owing to the divergence of 
the integrals, but we can apply it to the projected function C 

G R * G R l * C = C, (6.5) 
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under the only requirement that C^pojp) should satisfy assumptions (1) and (2) in the way 
a retarded function does and vanish rapidly enough at p — » oo to make the integral over 
G^oc (Po ) Coo (po ) convergent . 

Equation ( |6.4| ) is the equation of motion for Gr. In the X — > oo limit, it reduces to the 
well-known equation for Gr. For the Keldysh component of the propagator, the equation of 
motion is given by 

G R l * G K = Gr 1 * {hG R - hG A ) = 0, (6.6) 

where we have ignored the terms 0(e). Owing to the presence of the product G^ 1 * hG A , this 
equation cannot be verified directly (the integrals diverge). Instead (analogously to the case of 
the product Gr * G^ 1 ), one multiplies it from the left by the function C, which vanishes rapidly 
at p — > oo and satisfies assumptions (1) and (2) in the way an advanced function does. 

6.3 Resummed Schwinger-Dyson Series 

We write the Schwinger-Dyson equations in the form 

Gr = Gr + iGr *Yj R * Q r , Ga = Ga + iGa *Tj A * Qa, 

Qk = G K + iG R *T, k *Q a + %Gk *T, a *Q a + iG R *Y, r * Q k . (6.7) 

The formal solution (where all products are convolution products and the operators are kept 
in the proper order) is 

Qr = G R * (1 - iZ R * G R y\ g A = G A *{l-iY, A *G A )-\ (6.8) 



IK 



Q R * (% , u p )(G A l - G R l ) + iQ) * Q A . (6.9) 



To use the formal solution of the Schwinger-Dyson series, we assume that the functions 
G R (A), Qr{A), and T, R ( A ) satisfy requirements (1) and (2) for the retarded components in the 
upper and for the advanced in the lower semiplane. 

This assumption deserves a few comments: For the retarded (advanced) bare propagators, 
our assumption is valid. 

If the retarded component is real between the cuts on the part of the real axis, the Schwartz 
theorem tells us that assumptions (1) and (2) valid in the upper semiplane are also valid in the 
lower semiplane of the first Riemann sheet. 

At equilibrium, perturbation theory yields the full propagator as a set of Fourier coefficients. 
The analytic continuation in the energy plane is not unique. This freedom is used to choose an 
analytic continuation that satisfies requirements (1) and (2) defined in Sec. 3. The positivity 
property of the spectral density then implies that the propagator has neither zeroes nor poles 
off the real axistj. A further implication is that the exact self-energy S/j(po,p) at equilibrium 
also satisfies the properties (1) and (2). This is not guaranteed for approximate expressions for 
self-energy 

In the formal solution of the retarded propagator, the factors Gr and T,r alternate regularly. 
This fact can improve the convergence properties of some integrals. 
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Now it is easy to write down the resummed Schwinger- Dyson series for the retarded (ad- 
vanced) propagator (with any exact self-energy obtained by the perturbation expansion that 
satisfies our assumptions). In terms of the corresponding propagator calculated at X = oo: 

Gr(a),x (Po,P) = J dpoiPx {Po,Poi)GR(A),oo(Poi,P), (6.10) 

where 

Gr(A) >00 {Poi,P) = z. = } zrz^ - f zr. (6.11) 

1 - i^R(A),oo{P01,P)(*R(A),oo{P0l,P) 



Starting from Eq. ( |6. 1 1[ ) one obtains the HTL-resummedl^J' c3 retarded (advanced) compo- 



nent of the propagator without use of the gradient expansion. 

Some more work is necessary to calculate the Keldysh component. Now, in addition to 
individual terms, the sum Qr )00 (Pq,P) should also satisfy (1) and (2) i.e., the imaginary part of 
Sr,oo should be negative. 

However, there is a possibility that ImX# is positive in some kinematical region. Then 
the resummed Schwinger-Dyson equation for a retarded component can create the pole in the 
upper semiplane. However, this case is very questionable: one sums infinitely many retarded 
functions (i.e., the functions which vanish at t < t') and obtains the function which is not 



retarded (i.e., nonzero at t < £'). Such cases are usually classified as pathologye^ £2l). At this 
point one should cautiously consider the use of the "physical" gaugee3, in order to prevent 
eventual gauge artifacts. 

Some indication that, in some cases, Gr,oo(Po,P) does satisfy assumptions (1) and (2) comes 
from the HTL limit. Indeed, at equilibrium, the HTL limit ofGR,oo(Poi,P) must satisfy (1) and 
(2), as it is easy to verify. As the properties of density functions enter only through the thermal 
mass and the position of isolated poles, the same must be true of any distribution allowing the 
HTL approximation. 

Owing to the fact that the Keldysh component of self-energy does not satisfy the analyticity 
assumptions (1) and (2), we can only try to integrate expression ( |6.9| ) using approximate and 
numerical methods. 

However, it is possible that £# can be decomposed into two pieces satisfying assumptions 
(1) and (2) as retarded and as advanced functions, respectively: T, K = —1^k,r + ^k,a- F° r 



example, this happens in the case of one-loop self-energy Then Eq. ( |fc>.9| ) becomes 

Gk = Gr* {h{p , uj p )G^ + iE KiA ) *Ga~ iGr * {h{p , u P )G^ + iL K>R ) * Ga- (6.12) 

Owing to the fact that the functions Gr and Ga are not singular in the point po = zhuj p , the 
terms containing G^ 1 and G~^ cancel mutually. As one of the remaining convolutions includes 
factors of the same type (RR or AA), we are left with a single convolution multiplication. This 
convolution contains neither the advanced first factor nor the retarded second factor; thus, in 
general, it cannot be worked out in a simple way, and it will contain non-WTPF contributions. 
However, it may be performed at least numerically. 

The appearance of the non-WTPF contributions signals stepping out of the space of pro- 
jected functions. Indeed, the calculation of the more complex diagrams, containing subdiagrams 
resummed into Gk, will not enjoy advantages of the presented calculus. 
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Finally, we note here that the calculation starting with Eqs. ( |6.7|) and ending with ( 3.10J ) 



- ( p,12| ) cannot be performed with the true (i.e., calculated, in some miraculous way, to all 
orders) T, K , T, R , and S^. Indeed, we have anticipated that the true T> K , J2 R , and S^ contain 
non-WTPF terms, and thus one cannot use instead of G R and G^ their restrictions to the 
subspace of projected functions. Using the gradient expansion (under large- X assumption) 
one obtains familiar equations of motion for the Green functions of interacting fieldso B B ti3 . 
However, the advantage of the presented calculus will be observed in the properties of collision 
integrals, where one can expect considerable simplifications and the possibility of evaluating 
contributions of more complex diagrams. A more complete discussion of collision integrals is 
out of scope of the present paper, and we hope to publish it elsewhere. 

6.4 Pinching singularities 

The pinchlike contribution to the Keldysh component of the resummed propagator is expressed 
asc£l (we treat only the scalar case) 

G Kp = iG R * (-Ek,r + ^k,a) * G A , (6.13) 



where we have introduced the short notation T>k,r{A) = Mpo,^)^^) + ^k,r{A)- 

Similarly as in the case of resummed contributions, we can perform convolution between 
alike components (RR or AA). Then one can integrate the terms containing £# and E^r 
with respect to P02, and the terms containing S^ and Ej^a with respect to poi- The result is 
intriguing: 

1 .^ 1 



G Kp ,x (Po,P) = ~ dp 01 P Xo (p ,Poi)— 5—7 i^K- 



P01 ~^p + 2 ^epoi P01 -Up- 2i tPoi 

1 / „ e 2iX (p -poi) _ e -2iX (p -\w p ) _ ^ 

J2 A / dpoi —7- — r ^k,r{Poi + ie,P)- 



Ha^i V i7c(2p -p 01 - Xu p ) v ' ' ' p\ x - w$ + 2iep, 



01 



e -2iX (p -poi) _ e 2iX (p -Xujp) _ 1 \ 

+ / d Poi —77, x — s ^kXpoi ~ i^p)— — — . (6.14) 



The first term in Eq. ( |6.14| ) is a projected function (WTPF) that becomes the usual pinchlike 
term in the X — > 00 limit. It is this contribution to which the conclusions c3 about the 
cancellation of pinching singularities apply . However, the other terms are of non-WTPF 
nature; contrary to the case of the product of simple pole terms, the discontinuity along the 
real axis appearing in the functions Tir(a) and Y,k,r(a) now prevents the vanishing of these 
terms. 

A full discussion on pinching singularities in the finite-time-after-switching formulation re- 
quires more efforts and we hope to publish it elsewhere. 

7 Modifications of the Feynman rules 

In this section, in the framework of the generic field theory with bosons and fermions, we discuss 
the changes of Feynman rules that are due to the "finite time" assumption. We further analyze 
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the diagrams with respect to the question of energy nonconservation. Indeed, we find that this 
feature appears together with the non-WTPF contributions. 

The calculations performed so far already contain all of the modifications of the Feynman 
rules required by the finite-tj assumption. In coordinate space, the only modification is that 
the bare propagators [Eqs. ( f4.4|) and ( fQ>[) 1 are limited by < xq and < yo; thus they are 



projected functions. In energy-momentum space, the above change reflects in the change of 
propagators, vertices, and the overall factor. 

To transform to energy-momentum space, we choose some vertex j, arrange the orientation 
so that all lines i become outgoing, and use the propagators represented by Eqs. ( |4.4| ), (fOl), 
and (|4.6| ) (the pi momentum is joined to the line i). Exponentials attached to Xj are easily 
integrated with the help of Eq. (|5.9|) : 



2Wo dX]t ' 2n(-EtPi + ie)' [ ' 

After performing this integration, instead of the bare propagators we obtain their X — > oo 
limits [Eq. ( f4.6|) 1, which are the familiar propagators of the usual (£; — > — oo) theory. 

At the vertices the usual energy- conserving 5(J2iPoi) is substituted by i(27r) _1 (— J2iPoi + 
ze)" 1 . 

Under the momentum integrals there is a leftover factor at the vertices Ja (by subscript A 
we indicate that Ja are vertices with amputated legs): 

e -'E JA ^(E^^ A ^ A ) 5 (7 _ 2) 

where A = ± depends on whether the corresponding momentum is outgoing of or incoming to 
the vertex Ja, and ij A is running through the nonamputated lines. 

The overall factor in the case of two-point functions is treated in a simple way: introduce 
a slow Wigner variable as the average over the times of boundary vertices, and the relative 
time [Eq. ( p.l|) ]. Finally, one can Fourier transform over the relative time. There emerges 
an overall energy-smearing factor Px {po,Po) for two-point functions and similarly for n-point 
functions. In the case of n-point functions, the choice of variables is large and might not be 
unique; namely, depending on the diagram calculated, one chooses the most appropriate set of 
variables. 

The overall factor takes care of uncertainty relations: the larger the elapsed "time" Xq, the 
smaller the energy smearing. 

In the vertex factor the energy is not explicitly conserved. This energy nonconservation is, 
through the uncertainty relations, related to the finiteness of X . In the limit of infinite X , 
energy conservation is recovered. Here we want to argue that for some choice of propagators 
entering the vertex, the energy is conserved explicitly. To see this conservation, assume, for a 
moment, that at least one of the unspecified propagators (D, G, and S) related to the chosen 
vertex, say D^, is a retarded function, D R . In this case, one can integrate over q , close the 
integration path from above (owing to e lX o(Po,j+^) ^ dosing from below is out of question), and 
collect the contributions from singularities. If there are no singularities (and we know that 
conditions (1) and (2) are valid for bare propagators), one just obtains the energy-conservation 
condition S(J2i Pot)- The same is achieved with the outgoing momenta and advanced components 
of the propagator with closing the integration path from below. Now we are going to show that, 
indeed, one of these possibilities is realized at each vertex. 
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Each individual denominator (YliPoi ~ ^e) -1 (the lines are all oriented out) can be easily 
integrated. To demonstrate this, we have to sum over the indices of the corresponding vertex. 
We rename the basis (i,j),i,j = 1,2 into \p, v], where p, v — — 1 correspond to i,j = 2, and 
fj,,v = 1 correspond to i,j = 1. Then we find the relations of the type (we assume a three-field 
vertex, but the proof extends easily to any number of fields): 

D M = ±(D K -(iD R -uD A ). (7.3) 

The sum over the indices in the chosen vertex (S - , D -, G - propagators of the outgoing lines; 
the factor p for the negative coupling of the vertex to which the index-2 ends of the propagators 
are attached) is 

y^ l*>S\u,j\D\u,,p\G\u.A — jISrDrGr + (Sk + ^Sa)(D k + pDa)Gr 

+(S K + \S A )D R {G K + vG A ) + S R (D K + P D A )(G K + vG A )}. (7.4) 

Expression ( |7.4j ) contains only terms including at least one retarded propagator: Sr, or D R , or 
Gr. 

Thus one can integrate the terms separately and find that the factor (J2iPoi — ^e)" 1 is 
effectively replaced by i^S(J2iPoi)- 

As there is nothing special at this vertex (the indices \,p,v remain unspecified), one may 
conclude that this is a general feature. Nevertheless, one should do it very cautiously, step by 
step, while problems may appear at some degree of complexity. Then, as seen in Eqs. ( |5.10| ) and 



Q5.13Q , we find a new element in addition to the energy denominator (—poj +Poj+i — ze) _1 . One 
obtains the extra factor e~ lX °( poj ~ Poj+1 ~ lt \ With the help of this factor, even the contributions 
from the poles of the retarded component in the upper semiplane will decay exponentially with 
the time X . 

However, the diagrams with resummed self-energy subdiagrams are particularly sensitive. 
In this case, one is strongly advised to undertake an intermediate step: to Fourier transform 
the two-point function with respect to the relative time, to investigate the analytic structure, 
and then to perform the multiplication of two-point functions. Owing to the cuts in Qr< A ) and 
to the non-WTPF contributions to Qk, it is likely that this is just the point after which we 
have to live without the advantages of projected functions. 

8 Summary 

We consider out of equilibrium thermal field theories with switching on the interaction occurring 
at finite time (ti = 0). 

We study Wigner transforms (also in the relative time Sq) of two-point functions. To develop 
a calculation scheme based on first principles, we define a very useful concept of projected 
functions: a two-point function with the property that it is zero for x < U and for yo < tf, 
for ti < Xq and ti < yo, the function depends only on xq — yo- We find that many important 
functions are of this type: bare propagators, one-loop self-energies, resummed Schwinger- Dyson 
series with one-loop self-energies for the case of retarded and advanced components of the 
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propagator, etc. The properties of the Wigner transforms in the X — > oo limit are particularly 
simple if they satisfy these analyticity assumptions: (1) The function of po is analytic above 
the real axis (for a retarded component, but below it for an advanced component). (2) The 
function goes to zero as |po| approaches infinity in the upper (lower) semiplane. We find that 
these assumptions are very natural at low orders of the perturbation expansion. The convolution 
product of projected functions is remarkably simple, much simpler than what one would expect 
from the gradient expansion. 

The Schwinger-Dyson series, with bare propagators and self-energies being projected func- 
tions satisfying assumptions (1) and (2), is resummed in closed form without the need for the 
gradient expansion. The calculation of the resummed Keldysh component is simplified to a 
double (and under a certain analyticity assumption, to a single) convolution product. This 
contribution signals the stepping out of the comfortable space of projected functions. 

The Feynman diagram technique is reformulated: there is no explicit energy conservation 
at vertices, there is an overall energy-smearing factor taking care of the finite elapsed time (Xq) 
and the uncertainty relations. 

The relation between the amplitudes (valid at low orders of the perturbation expansion) 
of the theory with switching on the interaction in the remote past and the theory with finite 
switching-on time, enables one to rederive the results such as cancellation of pinching singu- 
larities, cancellation of collinear and infrared singularities, HTL resummation, etc. Previously, 
these results were considered applicable only to lowest-order contributions in the gradient ex- 
pansion. 

The question arises whether higher-order contributions also remain within the space of 
projected functions satisfying assumptions (1) and (2). The answer depends on the eventual 
positivity of Im£#, explicitly time- dependent perturbation, and the appearance of the one- 
loop approximated or resummed Keldysh component. The positive Im£# can create the pole 
in the upper semiplane in the resummed Schwinger-Dyson series. However, this case is very 
questionable: one sums infinitely many retarded functions (i.e., functions which vanish for 
t < t') and obtains the function which is not retarded (i.e., nonzero at t < t'). Such cases are 
usually classified as pathology t^' C3). The way of breaking the scheme explicitly is to introduce 
direct time dependence through time-dependent perturbation, or through the background field 
which depends directly on time: ifi(x,t) = <ft(t) + ip(x,t). In this way, one obtains the two- 
point functions which are not projected functions. A natural step out of the space of projected 
functions occurs in the calculation of the resummed Keldysh component of the propagator. The 
appearance of the non-WTPT contributions signals that the calculation of the more complex 
diagrams containing, subdiagrams resummed into Qk, will not enjoy advantages of the presented 
calculus. 
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Figure Captions 

Fig. 1: Finite switching-on time integration path. 
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